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We discuss the implementation of quantum logic in a system of strongly interacting particles. 
The implementation is qubitless since "logical qubits" don't correspond to any physical two-state 
subsystems. As an illustration, we present the results of simulations of the quantum controUed-NOT 
gate and Shor's algorithm for a chain of spin-1/2 particles with Heisenberg coupling. Our proposal 
extends the current theory of quantum information processing to include systems with permanent 
strong coupling between the two-state subsystems. 
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INTRODUCTION 

Current quantum computation and quantum informa- 
tion processing theory is based on the manipulation of 
two-state subsystems known as physical qubits 
(The physical qubits need not be identical to the binary 
digits used for encoding information, which are called log- 
ical qubits). One begins with a system of weakly inter- 
acting or noninteracting physical qubits and, in order to 
implement quantum logic gates, one either applies short 
pulses of an external field or "turns on" strong coupling 
between the qubits for short periods of time. For exam- 
ple, in a recent proposal Q, the logical qubit states |0) 
and |1) are represented by selected states of three spins 
in a chain. The states used have 7 = 1/2 and mi — 1/2, 
thus. 



|0>-;)5(in)-UT))IT), 



ID 



im>-;)g(in> + UT»iT>- 



(1) 



Quantum logic gates are implemented in this proposal by 
turning on the Heisenberg interaction between neighbor- 
ing spins for short time intervals. 

Physical systems with permanent strong interactions 
between the subsystems are not included in the current 
theory of quantum computation. We believe that there 
is a need to extend the theory to include these systems 
because such interactions can quite realistically be ex- 
pected in many of the types of physical systems that have 
been proposed as quantum computers. In this paper, we 
discuss the implementation of quantum information pro- 
cessing in these systems. As an example, we consider 
a spin-1/2 Heisenberg chain in a non- uniform magnetic 
field. We assume that the Heisenberg coupling cannot be 
turned off, and thus it is impossible to separate the sys- 
tem into weakly interacting subsystems. Nevertheless we 
show that qubitless quantum logic may be implemented 



using external electromagnetic pulses. Although this pro- 
cedure is not scalable to quantum computers with large 
numbers of qubits, these systems nevertheless can be of 
use in small quantum devices for opto-electronic inter- 
faces, quantum cryptography and communication, etc. 
These small devices would store quantum information 
and perform some elementary quantum manipulation. 



LOGICAL QUBITS WITHOUT PHYSICAL 
QUBITS 

We now discuss general requirements for the implemen- 
tation of quantum logic in a many-level system. We use 
each energy level to represent a number. For example, 
the numbers to 3 can be represented within a four-level 
system. Each stationary state of the system is identified 
with a binary number, 
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The subscript I denotes that this is a logical notation 
and the qubits here should not be identified with any 
physical entities. Since any quantum logic gate can be 
decomposed into a sequence of one-qubit rotations and 
two-qubit controUed-NOT (CN) gates |], we would like 
to be able to implement these gates on the logical qubits. 
For instance, a "7r/2" -rotation |Q on logical qubit "1", 
Ui = Ui(7r/2), corresponds to the following set of trans- 
formations, 
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and the modified controUed-NOT gate, CNio, would op- 
erate in the following manner, 

CNio|00)j = |00),, 
CNiolOl)/ = |01),, 
CNiollO)/ = 

CNioIll)/ = i\10)i. (4) 



Now, we assume that we can drive transitions between 
the energy levels by applying a cyclic potential, V, whose 
frequency is resonant to the energy difference between the 
levels. A transition between two states is possible only if 
the matrix element of V between those states is nonzero. 
We assume that any transition which "flips" one of the 
logical qubits is allowed. This assumption is valid for the 
model discussed below. If a complete transition between 
two states takes a time T to complete, then application 
of V for this period of time is known as a 7r-pulse. If a 
pulse is applied to one of the eigenstates for a time inter- 
val T/2 (7r/2-pulse), then a uniform superposition of two 
states will be created (assuming only resonant transitions 
occur). For example, upon application of a 7r/2-pulse to 
the state |00);, resonant to the energy difference between 
|00)/ and |10);, the transformation which results is. 



jth pulse, (with h = 1) is. 



|00)z^-^(|00)i+z|10)0, 
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and if the same pulse is applied to |10);, the transforma- 
tion is, 



|10),^i=(|10)i+z|00>;). 
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These are the first and third transformations in (|^). 
The second and fourth transformations can be imple- 
mented by the application of a similar pulse, resonant 
to the energy difference between jOl); and Thus, 
the one-qubit rotation, Ui, can be implemented by two 
resonant 7r/2-pulses. The CNio gate (Q), on the other 
hand, can be realized by the application of one 7r-pulse, 
resonant to the energy difference jlO); These 
are idealized implementations of the gates, as in reality 
non-resonant transitions occur, which may give rise to 
deviations from the ideal output. 



HEISENBERG CHAIN QUANTUM COMPUTER 

Consider a chain of L spins with uniform nearest- 
neighbor Heisenberg coupling, J. The system is sub- 
jected to a number of radio-frequency magnetic pulses. 
The Hamiltonian of the system under the influence of the 



H = H + V, 
if = -^c^fe4^-2J^Ifc.Ife+i, 
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where ujk is the Larmor frequency of the fcth spin, 1^ 
is the spin operator for the fcth spin, f2 is the Rabi fre- 
quency and Vj and <f>j are the frequency and phase of the 
pulse. Upon diagonalization of H, we obtain 2^ eigen- 
states. Except for two of these states (the ones with 
spin projection quantum numbers to/ = ±L/2), they are 
entangled superpositions of the individual spin product 
states. 

In the case of uniform Larmor frequency, both the to- 
tal spin and spin projection operators, and /f^j, com- 
mute with H . Thus the eigenstates have definite values 
of the total spin quantum number, /, and the spin projec- 
tion quantum number, m/. The effect of the application 
of a pulse, assuming only resonant transitions occur, is 
to give rise to transitions which change the spin projec- 
tion by Am/ = ±1, but leave the total spin unchanged. 
Starting from the ground state, which has / = i/2, only 
2/-I- 1 = L -I- 1 levels can be accessed using successive res- 
onant transitions resulting from the application of pulses, 
out of the total number of 2^ levels. 

Thus, in order to be able to access all possible lev- 
els, a non-uniform Larmor frequency (i.e., a non-uniform 
external magnetic field) is required. In this case, H com- 
mutes with I^^^ but not with l'^^^. Thus, the eigenstates 
have definite values of TO/ but one is not restricted to the 
subspace of states with a fixed total spin. 



CONTROLLED-NOT GATE 

As a simple example of the implementation of quantum 
logic in the Heisenberg system, we consider the modi- 
fied controUed-NOT gate, CNio, in a two-spin system. 
We label the eigenstates of H (expressed in terms of the 
eigenstates of /f and /|) with binary numbers, thus, 

|oo), = I TT), 

|01)/ = ai|Ti)+a2UT), 
|10>,--a2|n>+aiUT), 

|11>; = IU), (8) 

Where ai and a2 are functions of the ratio J /Soj, where 
5lj — LOi ^ uj{). In the limit J/Suj — > cxd, ai = a2 ~ 
l/\/2, and in the limit J/6uj — > 0, ai = 1 and 02 = 0. 
Note that in this scheme, not only are the logical states 
entangled superpositions of the spins, but the "0" and 
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"1" states of the logical qubits can't even be identified 
with any well-defined states of the overall system. Thus 
the scheme is qubitless in the sense that a logical qubit 
doesn't correspond to any sort of physical subsystem. 

The energies of the eigenstates are given by the follow- 
ing expressions: 



for the states |00)/ and jll); respectively, 
and for the states |01); and |10); respectively, 



(9) 



i^i,2 = +^T^. (10) 

The non-zero matrix elements of the pulse, V , in this 
system are, 

i{io\vmi = i{ii\vmi = -^e— *-^ (n) 

where f2o and fii are "effective" Rabi frequencies, 

^0,1 = (ai±a2)f^. (12) 

Therefore, all transitions involving a flip of one of the 
logical qubits are allowed. 

Using the Schrodinger equation, we can write equa- 
tions of motion for the probability amplitudes of the log- 
ical states, Coo, C'oij etc. In the interaction represen- 
tation, the evolution of each of the amplitudes will be 
coupled to the amplitudes of two other states. For exam- 
ple, the equation of motion for Cn is. 



^ ^ (13) 



If we apply a pulse resonant to the frequency of the tran- 
sition 1 10)/ ^ 1 11)/ 1 that \s V ~ ~ E2 then, ignoring 
the non-resonant transitions, the solution for Cn is, 

Cn(t) = Cn(0)cos + ^e-^Clo(0) sin 



(14) 



Similarly, the solution for Cio with this pulse is, 
Cw{t) = Cio(O) cos + le+'^CiiiQ) sin 



(15) 



Thus, a pulse with a duration 7r/f2i (which is known as 
a TT-pulse), causes a complete exchange of probabilities 
between |10)j and |11)/, while one with a duration 7r/2rii 
(7r/2-pulse) causes |10)/ to be transformed into a uniform 
superposition of the two states and likewise for |11)/. The 
TT-pulse of frequency v — E3 — E2 and phase (j> — is 
clearly what's required to implement the modified CNio- 
gate k). 



ERRORS IN THE IMPLEMENTATION OF THE 
CONTROLLED-NOT GATE 

In our preceding discussion, we assumed that pulses 
only drive resonant transitions. However, non-resonant 
transitions do occur with a finite probability, which can 
give rise to significant errors even in a single-pulse im- 
plementation of a quantum gate. Here we present the 
results of numerical simulations of the implementation of 
CNio. As an example, we begin with the following initial 
amplitudes for the logical states, 

Coo(0) = ^, Coi(0)-0, Cio(O)-^, and 

Cii(0) = 0, (16) 

and present the simulated final probabilities after a pulse 
of duration r = n/fli. The ideal output amplitudes are. 
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Coi(0) = 0, Cio(O) 
Cn(0) = ^. 



0, and 
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The simulated final probabilities are plotted in Fig. |i| 
as functions of the Rabi frequency for four different val- 
ues of 6lu. All parameters in the simulation are given in 
units of J. We chose loq = 100. One can see that the 
variation of the probabilities definitely shows a depen- 
dence on 6uj as well as an oscillatory dependence on the 
Rabi frequency. These dependencies can be explained by 
dividing the non-resonant transitions into two categories: 

• Near-Resonant Transitions. In this case there is 
only one such transition, namely |00); ^ |01)/ 
which, for Suj = 50, has a transition frequency of 
V = E\ — Eq = 99.98 as compared to v = E3 — E2 = 
98.98 for the resonant transition. 

• Off-Resonant Transitions. In this case there are 
two, |00)/ ^ |10)/ and |01)/ ^ |11)/ which, for 
Suj — 50, have transition frequencies of 151.02 and 
149.02 respectively. 

The effects of all of the above transitions can clearly be 
seen in the 6u! = 10 case. In Fig. 0(b), which shows the 
final values of |Coip, the large amplitude low frequency 
oscillations are the result of the near-resonant transitions 
and the small amplitude high frequency oscillations are 
the result of the off-resonant transitions. In Fig. |l|(c) for 
|CioP, the primary effect is that of the high frequency 
oscillations since this state is not directly affected by the 
near-resonant transition, rather it is affected by an off- 
resonant one. Figure |l|(d) for |Ciip again shows the 
high frequency oscillations. However, the oscillations are 
modulated and shifted by the near-resonant transition. 
This may be explained by a two-step process, jll); 
|01)i |00)i, an off-resonant transition followed by a 
near-resonant one. In Fig. Il|(a) for |CooP the probability 
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FIG. 1: Dependence of final CNio-gate probabilities, |Cfe;p on the Rabi frequency, Q, for different values of Suj. The values of 
Q and Slu are given in units of J. (a) |CooP- (b) |Coip. (c) |Ciop. (d) |Ciip. 



is primarily influenced by the near-resonant transition, 
but it is shifted upwards due to the off-resonant ones. 

In order to decrease the errors on the CNio, one can 
reduce the Rabi frequency, which reduces the effect of all 
non-resonant transitions. Another method of error re- 
duction can be explained by an analysis of Fig. ^ The 
graphs clearly show that increasing Slo decreases the in- 
fluence of the off-resonant transitions. This is due to the 
fact that it increases the difference between the frequen- 
cies of the off-resonant transitions on the one hand, and 
that of the resonant one on the other. When Suj is 250, 
the effect of the off-resonant transitions is almost negli- 
gible. 

The remaining error is due to the near-resonant tran- 
sition. This error can be reduced by making use of the 
"27rfc-method" ||, |j, which essentially involves choos- 
ing values of the Rabi frequency where the error due to 
the near-resonant transition is zero. That is, a Rabi fre- 
quency is chosen such that a 7r-pulse, while giving rise to 
the resonant transition will also cause an angular change 
of 27rfc (where k is an integer), due to the near-resonant 
one; so that the probabilities of the states affected by the 



near-resonant transition return to their initial values at 
the end of the pulse. The condition for this is, 

^Jnl + A2 T = 27rfc, (18) 

where r = n/fli as before, and A is the difference be- 
tween the near-resonant and resonant transition frequen- 
cies. Thus, for 6uj > 250, the errors can effectively be 
reduced by choosing values for the Rabi frequency satis- 
fying the 27rfc condition. (Here, we don't consider more 
sophisticated pulse shaping methods, which can also be 
used to reduce the errors in the implementation of quan- 
tum logic gates 0). 

SHOR'S ALGORITHM 

In order to illustrate the implementation of a quantum 
algorithm in our system, we have performed a simulation 
of Shor's algorithm for prime factorization We simu- 
lated a pulse sequence for the factorization of the num- 
ber four in a four-spin Heisenberg chain quantum com- 
puter. This is a 16 level-system and the logical states are 



5 



0.3 



0.25 



0.2 



o 



1^0.15 - 



0.1 



0.05 





0.005 



0.004 



CM 0.003 



O 



0.002 



0.001 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

n 



FIG. 2: (a) Probabilities of the logical states, \Cn\ , following the 41-pulse sequence for Shor's algorithm. The values of the 
parameters used in our simulations were, J = 30, too = 100, 5uj — 30, Q, — 0.5. The white bars indicate ideal output probabilities 
of Shor's algorithm, and the black bars are the results of our simulation, (b) The same probabilities, on an enlarged scale, to 
show the probabilities of unwanted logical states. 



10000)/, 10001);, 10010)/, ... 1 1111)/, or in decimal nota- 
tion 1 0)/, |1)/, |2)/, ... |15)/. We used the same sequence 
of quantum gates as the simulations on the Ising chain 
quantum computer ||^, but the pulses used were those 
appropriate to the Heiscnberg case. 

Shor's algorithm requires two registers of qubits. We 
designate the left two logical qubits as the x-register and 
the right as the y-register. We assume that the initial 
state is the ground state |0000)/ = |0)/. The algorithm 
then proceeds in three steps. The first step is to create a 
uniform superposition in the logical x-register. The state 
which results from this is, 
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(10000)/ + 10100)/ + 11000)/ + 11100)/) 
(|0)/ + |4)/ + |8)/ + |12)/), (19) 



of three 7r/2-pulses - for the transitions |0)/ ^ |4)/, then 
|0)/ ^ |8)/, then |4)/ ^ |12)/. The next step is to trans- 
form to the state, 



x=0 



(20) 



where y{x) — 3^ (mod 4). This is implemented using six 
TT-pulses. The final step is to perform a discrete Fourier 
transform on the logical x-register. This can be imple- 
mented using a sequence of so-called A- and B-gates, 
which in our system requires 32 (tt- and 7r/2-) pulses. 
The final state after implementation of Shor's algorithm 
ideally contains the following four states with equal prob- 
ability, 



The resonant pulse sequence to create this state consists 
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and 
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or, 

|3)/, |5),, and |7)/. (21) 

Figure ^ shows the resuhs of our numerical simulations 
of the 41-pulse sequence. The white bars correspond to 
the ideal values of the probabilities (1/4) of the states in 
(pT|). The black bars are the probabilities resulting from 
numerical simulations of the pulse sequence. The larger 
probabilities are shown in Fig. |^(a), and the smaller ones 
corresponding to errors are shown on an enlarged scale in 
Fig. |^(b) . It can be seen that the simulated probabilities 
for the states in (^) are very close to their ideal val- 
ues. The maximum deviation on these states is 0.0066. 
The probabilities of the unwanted states do not exceed 
0.0043. Four unwanted states have probabilities between 
0.0030 and 0.0043, and the remainder have probabilities 
smaller than 0.0005. The sum of the probabilities of the 
unwanted states is 0.016. 

CONCLUSIONS 

We have demonstrated a scheme for the implementa- 
tion of quantum information processing in a system of 
particles with strong permanent interactions. We envis- 
age that a number of these systems could be used for the 
design of quantum logical devices. We have also shown 
that a system of strongly coupled particles is qubitless 
- logical qubits don't directly correspond to any physi- 
cal two-state subsystem. Nevertheless, it is possible to 
implement quantum logic operations in this system. 

As an example, we presented results on the implemen- 
tation of the controUed-NOT gate and Shor's algorithm 



in a spin-1/2 Heisenberg chain in a nonuniform magnetic 
field. The results show that quantum logic can be imple- 
mented with resonant electromagnetic pulses and that 
the errors caused by nonresonant transitions can be ef- 
fectively reduced. 
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